THE COLLAPSAR AND SUPRANOVA MODELS 



CHARLES D. DERMER 

Code 7653, Naval Research Laboratory 
4555 Overlook Ave. SW, Washington, DC 20375-5352 USA 
E-mail: dermer@ gamma, nrl. navy, mil 

This rapporteur review summarizes results presented in Parallel Session GBT2 
(Gamma Ray Burst Theory 2) on the Collapsar and Supranova Models held 25 
July 2003 at the 10th Marcel Grossmann Meeting on General Relativity, Rio de 
Janeiro, Brazil. A central issue in GRB studies is the process whereby energy is 
released from the GRB engine. One scenario is the collapsar model, where the 
evolved stellar core promptly collapses to a black hole surrounded by a massive, 
intermittently-accreting torus of nuclear density material. A second scenario is 
the supranova model, where the first step of a two-step collapse process leaves 
behind a rapidly rotating neutron star stabilized by rotation, which later collapses 
to a black hole while making the GRB. In the supranova model, a powerful pulsar 
wind lasting days to weeks after the supernova makes distinctive signatures from 
the heating of the supernova remnant (SNR) shell, which should be discovered 
with Swift. This model also predicts nonstandard reddened excesses from the 
cooling SNR due to the range of delays between the two collapse events, contrary 
to observations of SN emissions from low redshift GRBs which favor the collapsar 
model. The observational basis for both models is critically reviewed. Problems 
of the collapsar/intcrnal shock scenario powered by the Blandford-Znajek process, 
and a two-collapse supranova scenario powered by the B-Z process or the Penrose 
mechanism, or even Ruffini's pair electromagnetic pulse, are briefly discussed. In 
view of the standard 7-ray energy upper limit, an external shock model involving 
jetted relativistic ejecta powered by an explosive event provides the most consistent 
explanation for GRB observations. An outline of a GRB model is proposed where 
the second collapse takes place within minutes to hours after the primary SN 
event, during which time loss of centrifugal support along the rotation axis of the 
neutron star provides a relatively baryon-clean polar environment along which a 
newly formed, rapidly spinning black hole drives collimated baryon-dilute outflows. 

1 Introduction 

The contribution^^ entitled "The Diagnostic Power of X-ray Emission Lines in 
GRBs," by Dr. Markus Bottcher (Ohio University, USA), and "Gamma- Ray Burst 
Progenitors confront Observations," by Dr. Davide Lazzati (IoA, Cambridge, UK), 
frame the current debate. Bottcher's article^ summarizes X-ray observations, and 
the case is made for a supernova (SN) taking place weeks to months before the 
GRB, so that the circumburst medium (CBM) contains dense ( £ 10 10 cm -3 ), 
distant (« 10 15 - 10 17 cm) SNR shell material whose illumination could account 
for the appearance of X-ray features in GRB prompt and afterglow radiation. These 
features give important clues about the GRB environment, and would be compelling 
if their significance were only greater. The marginal (3 — 5 er) detections near 
threshold make for concern about the reality of the features. Swift should clear this 
issue up, with its 24' x 24' FoV X-ray telescope reaching 2 x 10 -14 ergs cm~ 2 s~ 4 
in a 10 4 second observation, with autonomous slewing capability to get on target 
within one minute(!). 

At the time of this writing (April, 2004), it seems probable that all GRBs have 
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associated SNe as identified by optically thick photospheric SN emissions appearing 
as reddened optical excesses in GRB afterglow light curves some 10 - 30 days after 
the GRB event. In low redshift GRBs, and most clearly in the case ofGRB 03032£ 3 , 
the SN emissions are well fit by SN 1998bw light curves K-corrected and stretched 
to match the redshift of the GRB source. These emissions are powered by trapped 
shock energy and radioisotope heating (principally due to the 56 Ni^ 56 Co^ 56 Fe 
decay chain formed by explosive nucleosynthesis or accretion disk winds-). The 
optical follow-on observations indicate that any delays between collapses to a neu- 
tron star and black hole (or strange star) last for less than a few days. This places 
constraints on the angular mom entum and mass at formation, and the equations 
of state of GRB progenitors . 

Lazzatl^ argues that the issue of SN excesses in late (10 - 30 days) time optical 
afterglow light curves is settled, and he examines the low redshift (z = 0.17) GRB 
030329 in detail. As yet, SN excesses have not been observed from a GRB which 
also displays X-ray features. One possibility is that X-ray features will be detected 
at high significance by Swift from a GRB with a nearly coincident (within ,$ a few 
days) SN, in which case the original supranova model-'-', with delays between the 
two collapse events ranging from weeks to months, is ruled out. If Fe features are 
firmly detected, then this might imply ejection of material from an Fe core (which 
then does not explain what powers the SN Ic light curve). Another possibility is 
that the GRBs displaying Fe features form a disjoint set from those which display 
SN emissions. In this case, two separate population tracks for GRBs which produce 
SN emissions and those which produce X-ray line features from the illumination of 
a distant reprocessor are impliecffl It is also possible that the X-ray features are 
artifacts of analysis. 

As is apparent, the observational situation, particularly with regard to the X- 
ray emission lines and absorption features and edges in GRB spectra, is murky. 
Then what about the theoretical situation? 

The relativistic blast-wave external-shock model with both forward and reverse 
shock emissions is generally accepted as the basis to interpret afterglow observa- 
tions, and even to apply during the prompt phase in those GRBs which display 
enhanced optical reverse shock radiation. At least, this model provides a service- 
able tool to interpret radiation fluxes and to derive physical and environmental 
parameters (e e , es, jicbm, and to distinguish uniform vs. windfed CBM). The ra- 
diation processes that make 7-rays during the prompt (w £90) phase is an open 
question, with external and internal shock models each proposed. An external 
shock model, powered by a single explosive event, is favored on the basis of the 
data, as argued below. 

In this article, supranova-model predictions of signatures of SNR shells heated 
by plerionic emissions emitted by a newly born, rapidly rotating and highly magne- 
tized neutron star are presented. The heating of the SNR, and the interaction of the 
GRB emissions with this shell of material imply many predictions for the supranova 
model. We critique the collapsar and supranova models and argue that they are 
both inconsistent with the data (even leaving aside the X-ray observations). The 
observations do not (yet) show the predicted plerionic signatures of a supranova 
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model, nor do they imply standard reddened emissions at late times if a fraction 
of GRBs have delays of weeks to months between the two collapse events (that 
fraction, and the validity of the original supranova model, continues to decline in 
significance as more observations of reddened excesses accumulate). The collapsar 
model cannot explain the prompt GRB emissions, durations, efficiencies, nor the 
standard maximum-energy reservoir result. It predicts a wind- type CBM, which is 
only seen on rare occasion^. The baryon-dilute jet must be able to penetrate the 
massive stellar envelope without substantial mixing and reach T factors ^> 100. 

Some new ideas are needed. One possibility is that energy is dissipated by 
general relativistic effects within the ergosphere of a charged, rotating black hole 
as it is formed, rather than extracted when a magnetized disk is accreted onto 
a rotating black hole, as in the collapsar modeM Numerical simulations of core 
collapse of massive stars show that material can be evacuated along the rotation 
axis of the stellar core, leaving a relatively baryon-clean environmen Within 
a few minutes to 1 - 2 days of the SN, a second collapse would drive a jet in 
the poleward direction of the newly formed, rapidly spinning black hole to make 
the GRB. This model would make jet formation more feasible than in a collapsar 
scenario where nearly baryon-pure ejecta, whether in the form of a pair fireball or 
a Poynting-flux dominated flow, has to "drill a channel" through a massive stellar 
envelope and avoid the processes that are so effective at mixing 56 Co in SN 1987A, 
namely strong Rayleigh- Taylor mixing instabilities along the jet head and Kelvin- 
Hclmholtz instabilities along the jet spine and in the Rayleigh- Taylor fingers. 

Such a scenario does not depart from the well-established connection of long- 
duration GRBs with high-mass stars born in star-forming galaxies that end their 
lives as Type Ic SNe. Indeed, rapidly spinning Wolf-Rayet progenitors stars are the 
most likely GRB progenitors, as in the collapsaiP^ and some versions (though not 
of the supranova models. The main difference is in the nature of the collapse 
process and the origin and meaning of the prompt GRB emissions, with important 
implications for the physics of stellar evolution and black-hole formation. Whether 
the GRB is formed by an engine or explosion is important for understandin g th e 
physics of black hole formation. Here we argue for an explosive event (cf. Ref. 

Such a model brings into focus the "stand ard- energy reservoir result ,: ^^or, more 
accurately, the "standard-energy upper limit'^ISl and avoids many unpleasant issues 
that are often avoided in discussions of the collapsar scenario, such as low efficiency 
and absence of curvature signatures that are expected from colliding shells in the 
internal/external scenario, and the frequent presence of a rather uniform CBM as 
inferred from afterglow model fits to the external shock emissioni^l I also comment 
on current models for energy extraction from a rotating black hole, including the 
Blandford-Znajek effect and the Penrose mechanism. 

Both the coll apsar and supranova models have difficulties to explain the GRB 
phenomenon. W o ' ' have argued that the pulses in a GRB light curve are flashes 
made when a relativistic blast wave intercepts density inhomogeneities in the CBM. 
Delayed low-level radiations from fallback onto the newly formed black hole are 
possible (and this could help explain X-ray features should they turn out to be 
real). The GRB light curves are compatible with an explosion or a magnetar 
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dischare The relativistic ejecta shell can support a reverse shock structure to 
explain prompt optical and early radio afterglows. 

Reviews of GRB physics and observations that focus on observation:]^! an( j 
theorjEl both from the p erspective of the external shock modeP^and the inter- 
nal/external scenariJ22E3l can be consulted for further background. 

2 The GRB/Supernova Connection 

The discovery that the GRBs which are observed today took place at cosmological 
distances led to the development of the relativistic fireball/blast-wave model^that 
generalizes the theory of SNe to cosmic explosions with relativistic ejections. The 
coasting Lorentz factors T of G RB outfl ows reach values of 10 2 - 10 3 , as quantified 
by 77 attenuation calculationJ23232ZI Associated with this discovery are other 
important observational result that tie the long-duration (tgo ^, 1 s) GRBs to 
high mass stars and thence to SNe or some subset of SNe: 



GRB host galaxies have blue c olors, consistent with galaxy types that are 
undergoing active star formatio n 

GRB counterparts are found within the optical radii and central region s of the 
host galaxies, consistent with locations within star-forming spiral arm 

In some cases (e.g., GRB 011121), the surrounding medium is better fit with 
a wind-fed than uniform 

• Lack of optical counterparts in some GRBs may be due to high column den- 
sities and extreme reddening from large quantities of gas and dust in the host 
galaxjl^AI, as found in molecular cloud complexes where massive stars are born. 

• In other cases, however, the HETE-2 mission has shown that the "dark optical 
GRBs" are due to the large range of temporal indices of GRB optical light 
curves, which are sometimes very steep and rapidly fading^! 

• GRBs are associated with nearly contemporaneous Type Ic SNe in low-redshift 
bursts where detection of the SN emissions are possible (see Section 4); 

• GRB outflows are in the form of highly collimated jets. The evidence for 
beaming is inferred from achromatic beaming breaks in optical afterglow light 
curves that occur when the Doppler cone of the decelerating relativistic outflow 
is about the same size as the opening-angle of the j eiP. 

From the beaming breaks, Frail et al!^ inferred that the mean solid angle 
subtended by observed GRB jets is about l/500th of the full sky, so that the mean 
beaming factor /& = 0.2%. As a consequence, there are many hundreds of GRBs 
events that take place for every one that is detected. The absolute total ener gies 
in 7-rays from most GRBs cluster within a factor of ps 2 of ~ 1.3 x 10 51 erga^, 
typical of the kinetic energy in SN ejecta. By performing the statistics of GRB 
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sources, ond^H finds that the rate of GRBs in the Milky Way could reach w 10% of 
the rate of Type Ib/c SNe. 

The value of fb is currently a subject of debate and is in any case not precisely 
defined for structured jets. Larger values of fb 0.01 - 0.1 are sometimes claimed, 
which makes the individual GRB event rarer and more energetic, even if the local 
time- and space-averaged power into a star-forming galaxy remains about the same 
based on the hard X-ray/7-ray output, » 10 44 ergs Mpc~ 3 yr . Thus, not all Type 
Ic SNe harbor a GRB, and the fraction of Type Ic SNe that host GRB explosions 
may be as small as 1%. Nearby GRBs are likely to radiate less energy in 7-rays 
relative to ejecta kinetic energy compared to more distant GRBsPl. Rapidly fading 
GRBs also tend to be underluminous in 7-ray energjEl 

These lines of evidence suggest that GRBs are produced by collimated relativis- 
tic plasma ejected during an explosion or an active engine phase of a subset of Type 
Ic SNe formed from the collapse of the cores of high-mass stars. 



3 X-Ray Features and the Supranova Model 

This section critiques M. Bottcher's contribution^. See his article for more details, 
including additional references to the original literature. 

X-ray features have now been detected in at least 8 GRBs. Table 1 of Ref. ^ 
summarizes these observations, but the most important results are: 

1. low significance X-ray Fe K feature observed from GRB 970508^ in a sec- 
ondary outburst with Beppo-SAX, consistent with an 56 Fe Lya emission line 
blue-shifted by 0.1c from the redshift (z = 0.835) of the GRB source; 

2. low significance X-ray Fe K feature observed in GRB 970828221 ( z = 0.958) 
with ASCA, which is only consistent with the redshift of the GRB host galaxy 
if the feature is identified with an Fe K recombination edge in a nonequilibrium 
situation; 

3. X-ray emission features observed in the afterglow spectra of GRB 991216 (z = 
1.00) with ChandrJEI identified with H-like Fe and a recombination edge of 
Fe, along with weaker evidence for a recombination edge and emission line of 
H-like S, which originate in an outflowing moving at ~ 0.1c with respect to the 
source redshift; 

4. variable and declining Fe abso rpti on lasting for the first « 14 s of the 7-ray 
luminous phase of GRB 99070*^1 at z = 0.8435; 

5. an X-ray feature observed from GRB 00021 w ith Beppo-SAX, which had 
no optical transient counterpart or redshift measurement, but if identified with 
Fe Ka, corresponds to z = 0.47; 

6. multiple high-ionization emission features of S XVI , Si XIV, Ar XVIII and 
Ca XX (but no Fe Ka emission line or edge) observed from GRB 01121lESI 
(z = 2.14) with XMM-Newton (whose reality has been questioned); 
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7. low significance high- ionization Si XIV and S XVI blueshifted by 0.12c with 
respect to the host galaxy r edshift (z = 1.254) observed from GRB 02081^2 
with XMM-Newton from a GRB triggered by HETE-II; and 

8. low significance high-ionization S and Si (no Ni or Fe) emission features ob- 
served with XMM-Newton from GRB 030227 triggered by the INTEGRAL 
Burst Alert SysterrM The emission lines appeared in one segment of the 
observation, and the X-ray spectrum was lineless in an earlier segment. The 
redshift was not measured directly, but equal to z — 1.32 and 1.34 for Si XIV 
and S XVI, respectively, which when transformed from the emission to the 
observer frame implies z w 1.6. 

These results, taken individually and as a group, indicate that GRBs take place 
in environments highly enriched in metals, such as active star-forming regions. A 
highly metal-enriched and dense shell of material surrounds the sources of GRBs, 
and this shell is sometimes moving at speeds of rs 0.1c with respect to the source. 
A major concern is that although the significance is rather marginal, the amount of 
line energy is uncomfortably large and approximately constant, ranging from about 
6 x 10 48 ergs per line in GRB 030227 to » 2 x 10 49 ergs for GRB 970508. The line 
fluorescence is typically detected over observing time scales of « 10 4 s. The line 
energies in the sample of GRBs have tended to drop by a factor w 4 in the later 
bursts, but this trend may be misleading. Another oddity is that although 56 Fe 
Ka emission was reportedly observed in the first four GRBs, only emission lines 
associated with Si, S, Mg, and Ca were seen in the later GRBs when XMM-Newton 
was used, which may however reflect the different capabilities of the respective 
X-ray telescopes. 

Taking the observations at face value leads to some immediate requirements for 
models. 

If the total 56 Fe Ka line energy is 10 49 i?49 ergs, then the number of Solar masses 
in 56 Fe if each ion is individually illuminated once is 

M( 56 Fe) = 10 49 £ 49 ergs ^ 6E 49 ( 

M Q 5F 5 (1 M Q /56m p ) ■ 13.6(26 2 )Zf 6 eV F 5 Z| 6 ' [ ' 

where F = 5F$ is the fluorescence yield. That's a lot of Iron, which if made through 
stellar explosive nucleo synthesis, must already have decayed ( 56 Ni decays through 



electron capturd^SEl so the decay timescale could even be longer than w 100 
days). 

Type II SNe are thought to release perhaps w 0.07M Q of 56 Ni, and Type la as 
much as a Solar mass. Even if Type Ic SNe release « 0.5M Q of 56 Ni, multiple ion- 
ization episodes from a dense medium are required to account for the extraordinary 
energy and equivalent widths of the emission lines and absorption features. 

Another constrain1@SEEl applies to the size scale of the emitting region, which 
is assumed to be moving at most at mildly relativistic speeds. From the duration 
time of the line observations, given in units of 10 4 t4 s, then 

ri n \ r> < Celine 3 X 10 14 t4 . . 

(1 - COsd s hell)R ~ ~l+~Z TTz Cm ' ^' 
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where O s h e u is the angle between the fluorescing or absorbing shell medium and the 
observer. If the target material is moving in the transverse (equatorial) direction at 
the speed of 0.1i>_ic, a torus geometry is natural. This implies that the maximum 
duration between a first and second collapse event is ;$ v~\ day, assuming that 

the ionizing radiations travel near light speed. Vietri et al.ESl and Bottchen^l and 
FryenEO consider a model where an expanding shell materially impacts and shocks 
a torus of material. This geometry permits only short time delays between the shell 
ejection and the GRB and X-ray line production. The target shell material could 
be from an intense progenitor wind formed during the common envelope phase of 
a compact object and helium star. 

A distant reprocessor could be « 10 15 - 10 17 cm away and contain clumps of 
highly enriched material with densities ^ 10 10 cm~ 3 . Variable X-ray absorption is 
essentia lly u ndetectable from a GRB surrounded by a uniform CBM, even if highly 
enrichecPSl Extremely clumped eject cl^l or blueshifted resonance scattering in a 
high- velocity outflow^l can account for the variable absorption in GRB 990705. 
In this geometry, a shell of material is illuminated by the GRB emissions and 
could originate from a SN taking place months earlier. The target material would 
have locations and size scales corresponding to requirements of the external shock 
model for the prompt phase^l There are ranges of ionization parameters that are 
compatible with the observations for a geom etry-dominated model, though not an 
engine-dominated (i.e., collapsar) model^. 

Considerable heating accompanies multiple recombinations and reionizations, 
thereby reducing the recombination rate coefficient, so this model is criticized to 
make such a large line energy. Pair production could en han ce the line production 
efficiency. Other possibilities includes a funnelgeometrjE3 with a long-lived pho- 
toionizing source, which has not been detected^, and a pair screeri^to provide a 
local flux of UV/X radiation. The question of whether photoionization models or 
thermal models (or both) apply remains an open question. 

Taken together, the earlier ejection of a stellar envelope would be advantageous 
to explain the X-ray features, as in the supranova modeH This reasoning, given 
also the results in Chapter 2, leads to the conclusion that GRBs are a subset of SNe 
powered by a two-step collapse to a black hole of a massive ( ^ several M ) core 
stripped of its H and He envelope. The collapse event drives a subrelativistic outflow 
to power the SN, and a relativistic outflow to power the GRB. The GRB outflow is 
likely to be highly collimated. The excitation of torus material or the illumination 
of shell material in the polar direction could pro duc e the line signatures. Possible 
stell ar prog enitor candidates are Wolf-Rayet stara^J or merger events involving He 



core: 



!'»■»■■ »' 



If the X-ray identifications of Fe features are real, then GRB SNe must yield 
IM© of 56 Fe, which is a challenge for any model. 



4 Delayed Reddened Excesses and SN Light Curves 

... TO 

This section critiques and updates D. Lazzati's contribution , and summarizes 
evidence for SN emissions in the light curves of GRBs and implications for GRB 
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source models based on the examples of GRBs 980425 and 030329. 

GRB 980425, discovered with Beppo-SAX, was a smooth-profiled GRB with 
mean duration ~ 20 s, and is likely associated spatially and temporally with the 
Type Ic SN 1998bw in a spiral host galaxy at a distance of w 40 Mpc. Its apparent 
isotropic 7-ray energy was only w 10 ergs. It displayed very intense radio emis- 
sions, indicative of a related mildly relativistic outflow. GRB 980425 was probably 
an off-axis event with a much lower 7-ray to kinetic luminosity than in other GRBs. 
The kinetic energy is estimated from the radio luminosity using minimum energy 
arguments, or by multiplying the X-ray luminosity with time at a fixed observing 
time (e.g., 8 hours after the GRB) when the electrons are in the fast cooling regime 
and making a correction factor for the number of ions per electrode. Th e ki netic 
energy inferred from the radio luminosity of SN 1998bw was ^ 10 erg and 
exceeds 10 52 ergs when modeling the SN light curv These energies are similar 
to the kinetic energy of the Type Ic SN 2003L, but orders of magnitude greater 
than the radio luminosity in the Type Ic SN 2002ap. 

Comparison of GRB 980425 with the low redshift (z = 0.17) GRB 030329 
associated with SN 2003dh is instructive. Based on the jet-break times, typical 
bea ming-c orrected GRB energies have a "standard" 7-ray energy reservoir of ~ 10 51 
erg d ' ' (see Section 5). GRB 030329 displayed an early optical jet break at w 0.4 
d, and a much later 15 GHz jet break at « 15 d. In contrast, GRB 030329 has 
a 7-ray energy E 1 w 2 x 10 erg and is therefore underenergetic compared 
to a "typical" GRB (as also are XRF 020903 and GRB 980425). But in each of 
these cases, the kinetic energies E^ e ^ 10 51 ergJSal. (By contrast, GRB 031203, at 
z = 0.1, which is very similar to GRB 980425, has E 1 « 10 50 ergs and E ke w 2x 10 49 
ergs.) 

All this evidence points to a highly structured jet, or a "spine in sheath" ge- 
ometry (as in the case of 3C 273^) to account for GRB detectability and phe- 
nomenology. A dirty/clean standard top-hat jet fireball model is too simplified. 
Even if one wished to maintain a unifying approach to GRB observations, then 
angle-dependent effects from a structured jet must be considered. In any case, the 
number of GRBs hosted by Type Ib/c SNe should be < 10%. 

The discovery that a SN emerges in the afterglow of GRB 030329 tells us that at 
least a fraction of SNe host GRBs, and that SN and GRB explosions are most likely 
simultaneouJ^or, at most, separated by a few days. Here summarizes the evidence 
for delayed reddened excesses from SN emissions in GRB optical afterglows: 

1. GRB 970228 at z = 695, which reached a flux at v = 3.8 x 10 14 Hz at » 38 
d after the GRbPEH of F v = 0.87/iJy, corresponding to a stationary-frame 
luminosity of Aird\vF v {t = 25 d) = 10 43 ergs s" 1 (d L = 1.4 x 10 28 cm), about 
the same as the SN 1998bw template. Reicharll^ argues that this emission 
signature cannot be explained with a dust echo. 

2. GRB 980326 with unknown redshifP^. The peak of the light curve at 6558 
A(R-band is 7000 A, 1.8 eV, y = 4.3x 10 14 Hz) appears at ~ 20 days and has 
a peak flux density w 0.4/iJy, or a vL v power of « 10 43 ergs s _1 if z ~ 1. 

3. GRB 980425 (z = 0.0085, or d w 40 Mpc) associated with the Type Ic SN 
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With a peak absolute luminosity ss 1.6x 10 43 ergs s 1 , it requires 



~ 0.7 M Q of 56 Ni to power the light curved". 

4. GRB 980703 at z = 0.9661, consistent with a SN 1998bw-type excess at w 10 
- 20 days after the GRB, or no SlsPl. 

5. GRB 000418 at z = 1.119, which had R-band emission exceeding w 1 /xJy 
40 days after the GRB. If due to jetted afterglow emission, this flux would 
correspond to a vL v luminosity « 3 x 10 43 ergs s _1 that could conceal a SN 
excess^ 1 . 

6. GRB 000911 with a host galaxy redshift of z = 1.06 and optical emission 
with luminosity L ~ 1.5 x 10 43 ergs s _1 after as 1 week. The addition of a 
SN 1998bw-type excess improves the fit, and an alternative model of a shock 
interacting with WR progenitors involves several free parameters^! 

7. GRB 011121 at z = 0.362 with a distinct reddened excess appearing at w 20 
days after the GRB and reaching a vL v power of w 2 x 10 42 ergs s _1 at 7020A. 
This delayed emission is well fit by a SN1998bw template emission, a nd th e 
agreement improves if GRB 011121 exploded ~ 3 - 5 days after the Sr>P3Hl 

8. GRB 010921 at z = 0.451, which showed no SN to a limit 1.33 magnitudes 
fainter than SN 1998bw at the 99.7% confidence limilP^. 

9. GRB 020405 at z = 0.690, which showed a delayed reddened excess consistent 
with a SN 1998bw emissions dimmed by 0.5 mag nitude™ The radio data 



are not consistent with a density profile oc 



-2 



10. GRB 020410, redshift unknown, that exhibited an optical rebrightening one 
we ek aft er the GRB consistent with SN emissions from a nearly coincident 

snSE3. 

11. GRB 030329, asso ciated w ith SN 2003dh at z = 0.1685, with a highly variable 
173173131 



optical component 

12. The X-ray flash XRF 031203 at z = 0.1055, associated with SN 20031w, which 
can be represented by SN1998bw brightened by 0.55 maJSil! 

These observations show that in many cases, GRBs and XRFs are accompanied 
by a nearly contemporaneous SN which is spectrally much like the Type Ib/c SN 
1998bw that is temporally and spatially associated with GRB 980425. 

A collapsar model GRB must have a progenitor with stellar photospheric radius 
^ 10 11 cm, which limits allowed progenitors to massive stars that have lost their H 
and He envelopes, such as Wolf-Rayet stars or close binary systems^]. They must 
keep the ejecta extremely baryon-free by pushing away ~ 0.01Af Q while retaining 
a baryon purity ^ 10~ 3 Mq mi xed in with the explosion. A model for jet emission 
must explain the Amati relation^, where the uF v peak energy E p k oc , where 

Eiso,-/ is the apparent isotropic energy release. 
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Although Lazzati had argued that the bumps and wiggles in the light curve of 
GRB 02100^23 

were due to density enhancements in the CBM, he abandons that 
model for a model with episodes of delayed injection to explain the complex optical 
light curve of GRB 03032^23. He claims that 'any event that takes place over 
the entire fireball surface' will produce deviations in the afterglow light curve 'on 
timescales St ~ T, where T is the moment in which the deviation begins.' " This 
argument does not apply when the size scale of the density irregularities is much 
smaller than which is the most interesting regime from the point of view of 
an external shoc k mode l. Besides the complexity expected in the CBlvPl 

special 

relativistic effecti^"^ enhance emissions from shell-cloud interactions for clouds 

with e < i/r. 

The problem of the interpretation of the optical polarization curve of GRB 
030329 may be solved if complex behavior in its optical afterglow light curve is a 
result of localized blast-wave/cloud interactions in an external shock model. 



5 The Standard Energy Reservoir and the Collapsar Model 

Because of the potential importance of this finding, the "standard energy reservoir" 
result is summarized before enumerating the internal contradictions of the collapsar 
model. 

The standard 7-ray energy release from a GRB holds at least as an upper 
limit within a factor of ss 2, and was originally quoted as E 1 k, 3 x 10 51 /(?7 7 /0.2) 
er gjl4| with efficiency for 7-ray production compared to the total energy in 
the outflow given by r] 7 ~ 0.2. The absolute 7-ray energy from a GRB is given by 
E 1 — fbE ison , where the beaming factor h = 1 — cos9j, E ison — Andj^Fj/ (1 + z), 
and the jet-opening half-angle given by-^ 



9j = 0.057 



tj \ 3/8 / 1 + zy 3/s [ E ison 1 " 1/8 / „ 1 , 



\0.2J 



0.1 cm 



1 day J \ 2 J |_ 10 53 ergs 

(3) 

Note the weak dependencies on the unknown densities and efficiencies, which can be 
further constrained by afterglow modeling and comparison of the measured 7-ray 
energy with the kinetic energy of the flow. 

The standard 7-ray energy reservoir result is now superseded by observations^! 
of low E 7 outliers GRBs 980425, 030329, and 031203 which are under luminous, 
implying that the standard 7-ray energy is an upper limit and that there are biases 
favoring the detection of distant GRBs with maximum energy release. 

This result is quite difficult to explain in the context of a collapsar scenario, 
where the relativistic jet passes through a massive dense stellar envelope. A "large 
diversity in any accompanying SN components of GRBs is expected," and only a 
tiny fraction of a Solar mass of baryons mixed into the explosion would make for 
a considerably different amount of energy in 7-rays, s o th e "apparent constancy of 
the gamma-ray energy release is even more mysterious^! " 

Some other particulars which an internal shock/collapsar model must answer to 
be viable are: 
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1. How is it possible for pulses in GRB blast waves to have high efficiency in an 
internal shock model and also to have approximately constant (within a factor 
of 2 - 3) values of peak photon energy E p k in different pulses of the same GRB 
in spectra that extend to hundreds of MeV? 

2. If internal shells have such widely different efficiencies, then how is the standard 
7-ray energy reservoir result explained? 

3. Why are curvature effects from spherical GRB blast waves, a s foun d in an 
internal shock scenario, rarely if ever seen in GRB pulse proffieJSSHD? 

4. How is it possible for a jet to pass through a massive stellar envelope of several 
light seconds width and create pulses that are fractions of a second long? 

5. If the collapsar model is a failed SN, that is, there is no associated SN, then 
how is the enormous amount of 56 Ni produced to power the light curve^'/ 

6. Why have there been no detailed comparisons between temporal and spectral 
characteristics of GRB pulses, or even a smooth fast-rise, slow decay light 
curve, to show whether an internal shock scenario with high efficiency can 
explain GRB data? 

An external shock model with a single relativistic blast wave can simply answer 
these questions: 

1 . An external shock model is much more efficient than an internal shock mode^^, 
and approaches 100% efficiency for 7-ray production from accelerated electrons 
when the CBM is sufficiently dense. The values of E p k are relatively constant 
because the relative Lorentz factor in the interactions is about the same for 
interactions within the Doppler beaming cone. 

2. The most energetic GRBs are those where the largest fraction of the available 
kinetic energy is dissipated in the form of 7 rays. 

3. Curvature effects are not expected to be seen in an external shock model where 
the blast wave intercepts inhomogeneities on angular scales smaller than « 1 /Y. 

4. The jet originates from a bare collapse into a medium that has already been 
evacuated of much of the dense material. 

5. A two-step collapse process would first form a SN with a large amount of 56 Ni 
production, and would later collapse to a black hole to produce the GRB. 

6. Comparisons of pulse profiles, ti me his tories, and phenomenology with an ex- 
ternal shock model are 

favorablJSEEI. 



Zhang and Meszaroa=< y review the debate. 
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6 The Heated SNR Shell in the Supranova Model 



Any viable GRB model must be able to explain the increasing evidence for 
supernova-like emissions in the optical afterglows within the context of a model 
that explains the highly variable pulses in the prompt phase of a GRB. The evi- 
dence for X-ray features in prompt and afterglow GRB spectra indicates that high 
density, high mctallicity material near GRB sources. 

Is it possible to construct a model to explain at the same time the X-ray features 
and the delayed reddened excess emissions? Could a SNR shell that is illuminated 
and heated by a pulsar wind (PW) and cools following the GRB event make emis- 
sions that explains the delayed reddened excesses seen from GRBs? 

wJSni have tried to construct a model within the context of the two-step col- 
lapse process in the supranova scenario. This model has a number of attractive 
features: the period of activity of a highly magnetized neutron star preceding its 
collapse to a black hole can pr oduce a pulsar wind bubble consisting of low density, 
highly magnetized pair a"^" ' \ in accord with afterglow model fit*"4 The heating 
of the SNR shell by the X-ray synchrotron radiation from the pulsar wind electrons 
could provide an important source of ionizing radiation. A cooling SNR shell could 
produce the excess reddened emissions. 

We have found, however, that such a solution for the reddened excesses requires 
fine-tuning and, if not fine-tuned, predicts a wide range of cooling emissions that 
have not yet been observed. (It is important, however, to measure the spectra of 
SN emissions from nearby GRBs to determine how similar they are to SN 1998bw.) 
Thus the original supranova scenario appears not to be compatible with observa- 
tions of SN emissions if all GRBs, as seems probable, display rather uniform SN 
emissions in GRB optical afterglows. 

If the delayed reddened excesses in GRBs cannot be explained by the supranova 
model, it is ruled out. A simple derivation of the emissions from a SNR shell heated 
by a PW is presented in the Appendix. 

The results of that analysis are displayed in Fig. 1, which shows the total vL v 
flux composed of synchrotron, SSC, and thermal components. Standard parameters 
are used, as defined in the Appendix. The heavy solid and dashed curves correspond 
to standard parameters with t/t sc [ = 0.1 and 0.01 as labeled, where t s d is the 
magnetic dipole spin-down timescale, the crossing (or dynamical) timescale of the 
shell tdyn — 32.4 and 3.2 days, and T e fj — 1700 and 5420 K, respectively (t sc [ = 
3240 d). This system is in accord with the standard supranova model and might 
explain the delayed excess emission in GRB 011121, but requires fine-tuning. The 
set of three dotted curves employ standard parameters, except that B\2 = 4, and 
t/t sd = 1, 0.1, and 0.01 as labeled {t sd = 202 d), with t dyn = 20.2, 2.0, and 0.2 days 
(= 17 ks), and T eff = 3760, 11900, and 37630 K, respectively. 

The intense PW X-ray synchrotron emission provides a source of ionizing pho- 
tons throughout the SNR shell that would decay on the timescale tdyn (the line 
flux would decay on ss 2td yn due to emissions from the far side of the shell). The 
Bu = i,t/t s d = 0.01 case might correspond to a GRB 011211-type system. The 
PW radiation field provides an additional source of ionizing photons in the supra- 
nova model, and would be energetically important for this case if the typical en- 



main: submitted to World Scientific on February 2, 2008 12 




v(Hz) 



Figure 1. Total model spectral energy distribution composed of synchrotron, synchrotron self- 
Compton, and shell thermal emission from a pulsar wind and wind-heated SN shell. Standard 
parameters, given in the Appendix, are used, except where noted in the legends. 



ergies of electrons injected by the pulsar wind is 7 = 10 4 . A varying ionization 
flux impinging on an expanding SNR shell will produce a characteristic kinematic 
signature of shell illumination. 

The thin solid curve is the result for standard parameters except that B12 = 2, 
j w = 10 6 , and t = 0.lt sd = 81 d. For this case, t dyn = 8.1 d and T eff = 4360 
K. It is apparent that a range of parameter could make quasi-thermal emission 
between ~ 10 14 - 10 15 Hz at the level of ~ 10 42 - 10 43 ergs s _1 some tens of days 
after the GRB event, but only if accompanied by intense heating from the plerionic 
synchrotron radiation. 

Thus, a delayed reddened excess from the cooling shell may be seen, for plausible 
parameters, some tens of days after the GRB even if the GRB takes place months 
after the SN. The PW heating is in addition to any heating due to the Ni— ■> Co 
— > Fe chain, and may relax the need to have £ 0.5M Q of 56 Ni, as required to 
explain the unusually bright light curve of SN 1998bw^3. Whether the heated shell 
looks like the observations of reddened emissions in the optical emissions of GRBs 
will require better observations of Type Ic SNe and low-redshift GRBs, and better 
studies of this system. 

Should this model be correct, then GRBs would exhibit a bright hard X-ray 
emission feature corresponding to the plerionic component that could be detected 
with INTEGRAL or Swift from nearby z < 0.1 - 0.3 GRBs when the PW syn- 
chrotron flux radiates <L 10 44 ergs s _1 at X-ray energies. A nonthermal optical 
emission signature would accompany the reddened excess, with a relative intensity 
that depends sensitively on the energies of the PW electrons accelerated. A hard 
X-ray sky survey instrument such as EXIST could monitor for the PW nonthermal 
emission preceding aligned and off-axis GRBs. The X-ray synchrotron radiation 
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provides a source of external radiation to en hanc e photomeson production by en- 
ergetic hadrons over the standard GRB modepsl. 

The supranova model makes a number of other predictions: 

1. Compton-scattered 7 radiation fr om GR B electrons scattering photons from 
the intense plerionic radiation nelcEZIISJ 

2. High-energy neutrinos from the interactions of accelerated protons with the 
SNR shell, a nd due to photomeson production with the external radiation 

3. 7-ray emissions from the newly-formed, rapidly rotating pulsaJ^ prior to the 
second collapse to a black hole. 



7 Energy Release from a Newly Formed, Rapidly Rotating Black 
Hole 

At the heart of the GRB debate is the process that liberates an enormous amount 
of energy over a very short time in the form of a collimated plasma outflow with 
very small baryon contamination. Though some models have considered highly 
magnetized puls ar emissions^! or (as was discussed in the Session) collapse to 
quark matteJ^, it is generally assumed that GRBs are powered by the birth of 
a newly formed black hole. Until the question of the power source of GRBs is 
solved, problems about the origin of 56 Ni, which is supposed to power the reddened 
excesses observed in Type Ic SNc whether from explosive nucleosynthesis or disk 
windgSl — cannot be solved. 

Here we assume that a GRB is powered by the birth of a newly formed, rapidly 
spinning black hole. In addition to the enormous gravitational energy available 
during the collapse process, a spinning bl ack hole can be formed with ^> 10 52 ergs 
of rotational energy. Ruffini and colleagues^! argue that polarization of the vacuum 
can produce an e + -e~ pair-electromagnetic pulse of enormous energy that, given an 
external shock model, could explain the light curves of GRBs. The problem here is 
that the pair pulse is isotropic, contrary to inferences from observations that GRB 
ejecta are highly collimated. Moreover, unless the collapse to a black hole takes 
place well after the SN, it would be difficult for an isotropic explosion to break 
through the stellar envelope. 

Tapping the spin energy of the newly formed black hole seems a better bet for a 
GRB model. Advocates of a collapsar model generally invoke the Blandford-Znajek 
process where an induced voltage, formed when the intense accretion-disk magnetic 
field is shorted out in the ergosphere of the spinning black hole, accelerates particles. 
For a sufficiently massive and highly magnetized accretion disk surrounding a black 
hole with spin parameter a — Jc/GM 2 — > 1, adequate power can be injected 
to explain a Though this is the most popular model for GRBs at the 

present time, quantitative explanations of temporal and spectral features of GRB- 
like pulses, the constant energy-upper limit, 56 Ni production, the escape from the 
stellar envelope of a lage T outflow, the mean duration of GRB activity, etc. (see 
Section 5), remain unsolved. Another problem is that simulations of collapse to a 
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black hole do not leave behind a sufficiently massive torus to power a 

GRE 5nnsi 

Such simulations do not give much support for a supranova scenario, either, because 
differential rotation induces toroidal Alfven waves that can magnetically brake and 
destabilize the neutron star, thereby causing collapse within hundreds of seconds. 

Another possibility is that the Penrose process, invo lving either Penrose Comp- 
ton scattering or Penrose pair production ! ! 'M , extracts energy from the 
newly formed, rapidly rotating black hole. Particles on retograde orbits in the 
ergosphere can extract rotational energy through decay or scattering. It is not 
clear, however, that retrograde scatterers can be formed in the dynamically form- 
ing black-hole, nor is it clear that the efficiencies of this process are adequate to 
power a GRB. 

Given the uncertainties in our knowledge of the fundamental GRB energy release 
process ( "a black hole within a black box" ) , it might be best to speculate on a viable 
GRB model from the observations, as summarized earlier. 

1. The delayed reddened excesses are indeed due to SNe, and the progenitor star 
of a GRB collapses to a neutron star to form a Type Ic SNe, forming 56 Ni 
through explosive nucleosynthesis. 

2. The delays between collapses from a neutron star to a black hole last from 
minutes to days rather than from weeks to months. This means that X-ray 
features associated with Fe are mistakenly identified, though the signatures of 
S, Si, and Ca could still be real. 

3. During the period between the two collapse events, powerful pulsar winds 
could evacuate a channel along the rotational axis of the progenitor star and 
the pulsar, where the overlying material has the smallest column density. 

4. A single, highly collimated explosive event during the collapse of the neutron 
star to a rapidly rotating black hole could eject a relativistic outflow to produce 
the GRB prompt and afterglow emissions through external shocks with the 
CBM. 

Swift will be crucial to test this and other GRB models by providing better 
measurements of the X-ray features, by discovering many nearby GRBs from which 
delayed reddened emissions can be precisely measured, and by searching for signa- 
tures of a heated SNR shell. 

Acknowledgments 

The work of CD is supported by the Office of Naval Research and GLAST Science 
Investigation Grant No. DPR-S-1563-Y. 

Appendix 

A Pulsar Wind Physics 

A pulsar wind (PW) and pulsar wind bubble (PWB) consisting of a quasi-uniform, 
low density, highly magnetized pair-enriched medium within the SNR shell is formed 
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by a highly magnetized neutron star during the period of activity preceding its 
collapse to a black hole-- 5 * 1 . The interaction of the PW with the shell material will 
fragment and accelerate the SNR shell, and the PW emission wil l be a source of 
ambient radiation that can be Comptonized to gamma-ray energieJ^ZEl 

The nonthermal leptonic PW radiation provides a source of hard ionizing radi- 
ation to heat and photoionize material in the SNR shell, which can be described 
in terms of a partial covering model as a result of the strong PW/SNR interac- 
tion. This radiation heats the shell, and the escaping thermalized radiation is seen 
as the delayed reddened optical excesses in GRB optical afterglow light curves. 
The source of ionizing nonthermal radiation, which is extinguished when the GRB 
occurs, would provide a nonthermal precursor to a GRB, though at low flux. 

A very rapidly rotating, strongly magnetized neutron star is formed in the first 
step of the supranova scenario. For a neutron star with surface polar field of 
10 12 i?i2 G, a rotation rate of 10 4 fi4 rad s _1 , and a radius of 15ris km, the radiated 
power L = 7 x 10 44 P>l 2 r%nj ergs s _1 . The quantity E rot = \j{GM 2 /c 2 )n ^ 
2 x 10 53 jo.5-M|f24 ergs = 10 53 E^ is the initial rotational energy of the neutron 
star, expressed in terms of the neutron-star mass M — SM^Mq and the dimension- 
less angular momentum j'0.5 = GM 2 /2c, which remains roughly constant during 
collapseP. The spin-down time t sd = (E rot /L ) = 2.8 x 10 8 J0.5-M3 (.B? 2 rf 5 fif ) _1 s. 

The torque equation, (l = —Kfl n , where the braking index n = fltt/Ct 2 , can be 
solved to give fl(t) = fi (l + tAd)" 1/(n-1) , provided that K remains constant!^. 
(For the Crab, n — 2.5.) We assume that the spin-down power is radiated as a 
PW in the form of particle and Poynting flux, so that the wind luminosity L w = 
L s d = (£e+£,p + £,B)L s d = L Q /(l + t/t sd ) k , where k = (n + l)/(n-l), and we have 
divided the wind power into leptonic ("e"), hadronic ("p"), and electromagnetic 
("B") fractions £j, i = e,p,B. 

In the supranova model, rotational stability is lost due to spindown, and the 
neutron star collapses to a black hole to produce a GRB at time t co u. Here we 
assume t co u £ t s d, so that for k > 1 the wind energy E w (t) — r\ (k — 1) Lot s d [1 — 
(1 + {t/t s d) 1 ~ k } — > r/Lot when t ^ t s d, where the parameter 77 roughly accounts for 
the escape of wind energy from the SNR shell as well as losses of wind energy due 
to radiation, and is assumed here to be constant with time. The equation for the 
SNR shell dynamics in a spherical approximation for the SNR shell is M$nrR = 
AiTR 2 p w {t), where M$nr = itisnrMq is the SNR shell mass, the wind pressure 
p w (t) ~ E w (t)/W pwb , the PWB volume V pwb = 47ri? 3 (i)/3, and the PWB radius 
R(t) = J* dt' v(t'). This equation can be solved^ to give 



R{t) = v t(l + J - — ) , where t acc = S 1.9 x 10 7 s (4) 

V z acc IT/J^O r l n 12 r lb l ' l 'A 

It follows that t acc /t sd = Eg NR /r,E rot ~ O.lEf^/r?^*, where \M SNR vl = 9 x 



l0 51 m S NRP-i = 10 52 E£ 2 NR ergs is the SNR shell kinetic energy, and v a = 0.1/3_ x c 



J ke I'l^rot — "--^^52 / '/^53 ) ™ n,:lc 2 J 
T 

J 52 

is the SNR shell coasting speed. If rj = 1, then t acc -C t s d for standard values, 
and shell acceleration must be considered. A highly porous shell has 77 <C 1. This 
will occur for strong clumping of the ejecta, which takes place on the Rayleigh- 

Taylor timescald^ t^.T — yx/R, where the characteristic clumping size scale is 
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x = fR. Hence (tRT/tacc) = f 1 ^ 2 (t/tacc) 3 ^ 4 in the regime t £ t acc . This shows 
that small-scale (/ -C R) clumping will occur when / -C \/t/t acc . Only a detailed 
hydrodynamic simulation can characterize the porosity of the shell due to effects of 
the puls ar w ind, which 77 parameterizes. Here we assume, in contrast to the picture 
of Ref. p^j, that a large fraction of the wind energy escapes the SNR shell and 
77 ~ 0.1, so that shell acceleration can be neglected. This effect also causes the shell 
to become effectively Thomson thin much earlier than estimated on the basis of 
a uniform shell approximation, except for high density clouds with small covering 
factor. A general treatment is needed to treat shell dynamics, including travel-time 
delays of the pulsar wind and deceleration of the SNR shell by the surrounding 
medium. When t co u £ t acc , t s d, as assumed here, R = v$t. 

The strong PW provides a source of nonthermal leptons, hadrons, and elec- 
tromagnetic field. Dominant radiation components considered here are leptonic 
synchrotron and SSC radiation, and thermal emission from the inner surface of the 
SNR shell which is heated by the nonthermal radiation. 



A.l Wind Synchrotron Radiation 

The volume-averaged mean magnetic field B w in the SNR cavity powered by the P W 
is obtained by relating the wind magnetic field energy density ub — B^(t)/8ir = 
^BwU w (t) = ZeBw'nLQt/A'KR? through the magnetic-field parameter es w , also as- 
sumed to be constant in time. Thus 



B w (t) = - , where B = J § = 4 x lO*y/e B wV Sr^) G - s . (5) 

t V v ° P-x 

Let j w — 10 5 75 represent the typical Lorentz factors of leptons in the wind zone 
of the pulsar. This value should remain roughly constant when t ^ t s d- A quasi- 
monoenergetic proton/ion wind may also be formed, though we only consider lep- 
tonic processes here, and furthermore do not treat nonthermal particle acceleration 
at the PW/SNR shell boundary shock. When synchrotron losses dominate, the 
mean Lorentz factor of a distribution of leptons with random pitch angles evolves 
in response to a randomly oriented magnetic field of mean strength B according to 
the expression — jsyn — (<?tB 2 /6irm e c)~f 2 , giving 

1 W = [~-T B (±- h- 1 , (6) 
7 U t 

where U is the injection time, Tj = j w is the injection Lorentz factor, and 



r /; ,„ = -f^l 7 ^2x 10 13 < ,,.,,/ ( ^gfc^ ) 7.-. s . 171 
D7rm e c p_ 1 



Writing the nonthermal electron injection function as 
dN e (ji,ti) r,( e L 

the electron Lorentz factor distribution dN e ('y;t)/d'y can be solved to give 

7^^^(1 + 1-1)-. (9) 
«7 t 7 
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In this expression, M e = tj^Lq/ ' (m e c 2 ^ w ), and we introduce the dimensionless 
quantities t — t/Tgj w and 7 = j/"f w . Adiabatic losses, which are significant on 
time scales ~ i/3, are negligible in comparison with synchrotron losses of wind 
electrons when i<l, and we restrict ourselves to this regime. 

The vL v synchrotron radiation flux vh a J ln ~ ^UBCaT^N e {^/ s ) 1 where 7 S = 
yJe/eB, e = hv/m e c 2 , £3 = B/B cr , and the critical magnetic field B cr = 4.41 x 10 13 
G. Thus 



vLT{e) = ca T (At b J 2 J (j f' 2 * & + M - 1)~ 2 , for c < w , 

(10) 

where u = B/ lw B cr T B S 4.5 x 10- 19 /3 3 _ / 1 2 /(V^i^Bi 2 r 3 5 ^i7 5 ), and 

M 2 __^7i 4.4 x IP- 9 /3 3 / 2 7 5 
e _ _ £B7tu _ ^ 4 _ { ^__ B ^ nl . (11) 

The synchrotron cooling timescale for wind electrons is t syn — j w TBt 2 . In 
the regime t 0$ where the wind electrons strongly cool, we can approximate the 
synchrotron spectrum by 

vLT[e) s K eLo ( _^)i/a ff(e . £0j ^ 5 (12) 



where H(x; a, 6) = 1 for a < x < 5, and H — otherwise. The value of £ at t s j is 



Jf^_ s 1.4 x io-5 jgjjg^i 
Ib7u, e Bw r?Sf 2 r}|n^ 1 75 



Strong cooling holds when i < t sc i/TB"/ w , noting that time in physical units can be 
found from equation Q. 

A. 2 Wind Synchrotron Self-Compton Radiation 

The relative importance of synchrotron self-Compton (SSC) to synchrotron cooling 
is given by the ratio p = u S y n /u B of the synchrotron radiation energy density to 
Ub = B 2 /8irt 2 (see cq.jS]). The PW synchrotron radiation energy density u syn ~ 
n{vL™ ax ) j '4ttR 2 c, where vL™ ax = (, e Lo and the parameter k > 1 is a reflection 
factor such that k = 1 corresponds to perfect absorption or direct escape (small 
covering factor) of the PW synchrotron radiation, and k> 1 corresponds to a highly 
reflecting SNR shell with large covering factor. One obtains p — 0.033K^ e P-i/eBwil- 
For large porosity with k = 1 and r\ sa 0.1, we see that the SSC contribution is 
small compared to the synchrotron flux when £ e /3_i/ 'esw ~ 1- We restrict ourselves 
to this regime where equation © is valid. 

In the 5-function approximati on fo r the energy gained by a photon upon being 
scattered by relativistic electronfllHI with scattering restricted to the Thomson 
regime, the vL u SSC spectrum from the cooling wind electrons is given by 

-SSC ^ VTKteLp N e j £ r 1 „-2 w x xiixi 



8ttv 2 T b t 2 V e max [ a (a + x) K a + x nlxo1 y ' 
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where a = i 1 - 1, x = TwV^oA' and x l = 7™ v /max ( e ' e_1 > £max)/e, where we 
have used approximation H12|) for the synchrotron spectrum. 

A. 3 Thermal Emission from the Interior of the Shell 

The SNR shell is pictured as having been shredded by the PW and therefore highly 
porous, though permeated with s mall scale (~ 10 12 -10 14 cm) density irregularities 
consisting largely of metalJ^^I Note that O.lm56.-iM of 56 Ni provides a ra- 
dioactive decay power of ~ 2 x 10 42 tos6,-i ergs s" 1 before decaying on an w 120 day 
timescale. The partial covering factor p c (< 1) corresponds to the covering fraction 
by optically-thick, dense shell inhomogeneities. Recognizing the severe limitations 
in the following expression, we determine the effective temperature T e ff in the 
interior of the shell through the relation 

4nR 2 a SB T^ ff p c £ A Pc + 2 x 10 42 m 56 ,-i exp(-t/120 d) ergs s" 1 , (15) 

where A{< 1) is an absorption coefficient, dependent in general on the evolv- 
ing radiation spectrum, though A is here assumed constant. The wind heat- 
ing is more important than heating from the decay of radioactive Ni when 
(X/0.1)(p c /0.1) (&/1/3) £ TO 56 ,-i/-B 2 2 r? 5 ft 4 . In this regime, 

letting k w 1 because p c ~ 0.1. For standard parameters, the shell is heated to 
~ 10 4 K temperatures when t/t s d <C 1. When Byi ~ 4, this occurs at t/t s d ^ 1. 
The PW heat source is extinguished once the GRB occurs, and the thermal emis- 
sion from the SNR shell will decay on the dynamical timescale td yn — R(t)/c = 
324/3_i (jo.siVfg /S 2 2 rf 5 f2|) (t/t s d) d, noting that the cooling timescale of a hot 
shell is generally shorter than tdyn- 

The parameters used in Fig. 1 are jo. 5 = -B12 = ^15 = ^4 = = K = M3 = 
"152,-1 = 75 = I, £e = e Slu = 1/3, and p c = i] = A = 0.1. 
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